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Quantum gravitational Kasner transitions in Bianchi-I spacetime 
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Department of Physics and Astronomy, Louisiana State University, Baton Rouge, 70803 

Due to non-perturbative quantum gravitational effects, the classical big bang singularity is 
replaced by a quantum big bounce of the mean scale factor in loop quantization of Bianchi-I 
spacetime. An important issue is to understand various differences in the physical properties 
of the spacetime across the bounce. We investigate this issue in the context of various 
geometrical structures, identified by the Kasner exponents of the metric, which arise on 
approach to the singularity in the classical theory. Using effective spacetime description of 
Bianchi-I model in loop quantum cosmology with dust, radiation and stiff matter, we find 
that as in the classical theory, geometrical structures such as a cigar or a pancake form, 
but they are finite and non-singular. Depending on the initial conditions of the matter and 
anisotropies, different geometric structures are possible in the pre- and post-bounce phases 
^ . in physical evolution. Thus, quantum gravitational effects can cause a Kasner transition in 

Bianchi-I spacetime, which is not possible at the classical level. Interestingly, we find that not 
CO i all transitions are allowed at the level of effective dynamics in loop quantum cosmology. We 

find the selection rules and underlying conditions for all allowed and forbidden transitions. 
' ^j ' , Selection rules suggest that for a given set of initial conditions on anisotropies, occurrence of 

U" 1 Kasner transitions follows a distinct pattern, and certain transitions are more favored than 

| , others. 

<N 

> ■ I. INTRODUCTION 

CO 

, Singularities such as the big bang are the boundaries of spacetime in Einstein's theory of general 

relativity (GR). Signalling the limits of validity of GR, and generally characterized by inextendibil- 
ity of geodesies and divergence of curvature scalars, they come in a diverse variety. They can be 
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strong or weak (harmless) [H, 0], and can arise by divergence of Ricci or Weyl curvature, or both. 
They may have various geometries such as a barrel, cigar, pancake or a point, as in a Bianchi-I 
spacetime 0-0] • These geometrical structures are described via Kasner exponents in the space- 
time metric, and it has been conjectured by Belinski-Khalatnikov-Lifshitz (BKL), that a general 
approach to singularity may involve chaotic oscillations between different Kasner phases j3, @]- 
However, for the matter satisfying weak energy condition, such transitions require a spacetime 
to have both spatial curvature and anisotropies, hence, these are absent in a classical Bianchi-I 
spacetime. For such spacetimes, the geometric structure of the singularity is determined from the 
initial conditions on the matter content and anisotropies, through the resulting interplay between 
the Ricci and Weyl parts of the spacetime curvature. As an example, for a Bianchi-I universe in 
GR, filled with dust, only cigar and pancake structures are possible, whereas if filled with stiff 
matter, barrel, cigar and point like singularities can form [5(. 

Such singularities are expected to be resolved in a quantum theory of gravity. Though the exact 
mechanism of this resolution is still to be understood for arbitrary spacetimes in quantum gravity, 
analysis of various isotropic and anisotropic models based on loop quantum gravity (LQG) suggests 
that the singularities are replaced by a non-singular bounce which occurs due to underlying discrete 
quantum geometry, and bridges two classically disjoint branches of continuum spacetime. The 
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occurrence of bounce, first obtained for spatially flat isotropic models in loop quantum cosmology 
(LQC) for massless scalar field (which is equivalent to stiff matter in terms of equation of state) 
9^1lK has been established in various matter models 1214181] and in the presence of anisotropies 
19- 241. 1 For an up to date review of recent results, see Ref. [i^]. A robust feature of these 



models is that non-perturbative quantum gravitational effects lead to singularity resolution when 
spacetime curvature becomes Planckian, irrespective of the initial conditions. These effects occur in 
a very short range of spacetime curvature, and quickly die when spacetime curvature becomes small, 
roughly at a percent of the Planck value below which LQC approximates GR. Various physical 
insights on quan tum geometry have been obtained using effective continuum spacetime description 
of LQC 3CH321] . which provides an accurate representation of the underlying discrete quantum 
dynamics for states which describe a macroscopic universe at late times U, 16, 18, 23|, 33|, 3|. Using 
effective dynamics, in context of the details of singularity resolution, bounds on the energy density, 
expansion and shear scalars in different models have been found 35j , and strong singularities have 



been shown to be generically resolved in isotropic [36J, [37( and Bianchi-I models [38|]. For more 
discussion on semiclassical limit and effective dynamics see Sec.-V of Ref. [29J] and the references 
therein. 

These developments set a stage of posing various interesting questions which fall beyond the 
realm of the classical theory. Of these perhaps one of the most important one concerns with 
understanding the way physics in the pre-bounce branch influences the physics in the post-bounce 
branch (and vice- versa). For the isotropic models in LQC, useful insights have been gained by 
considering properties of semi-classical states when evolved across the bounce 3914411] . As the 
bounce is approached, the geometrical structure of the spatial geometry is always point like on 
both sides of the bounce in the isotropic and homogeneous models. These results can change 
in anisotropic setting where due to the presence of Weyl curvature, physics is much richer. In 
particular, given that in Bianchi-I spacetime, the classical theory leads to different geometric 
structures as the singularity is approached, it becomes important to understand how such structures 
arise in effective spacetime description of LQC and the way they are affected across the bounce. 
It is important to note that far away from the Planck regime, since there is little distinction 
between LQC and GR, for a given set of initial conditions on the matter and anisotropies, a co- 
moving observer in GR and another in LQC measure practically the same Kasner exponents of 
the spacetime metric. Such observers thus deduce the formation of the same geometric structure, 
such as a cigar or a pancake, in high curvature regime. The important difference between their 
description is, that the observer in GR, is inevitably annihilated as the geometric structure becomes 
singular at the big bang, whereas the observer in LQC finds geometric structure to be finite and 
singularity free. 

These and related aspects of new physics have so far remained unexplored in the studies on 
Bianchi models in LQC. 2 Our goal in this work is to investigate these issues for Bianchi-I spacetimes 
where detailed results on the formation of geometric structures at the classical level are already 
available 0]. We consider effective dynamics of Bianchi-I model in LQC for a barotropic fluid 3 with 
equation of state corresponding to stiff matter, dust and radiation, with an objective to understand 
answers to the following questions: Given a set of initial conditions under what conditions do 
particular geometric structures arise in LQC as one approaches the bounce of the mean scale 
factor? Does a physical evolution across the bounce yield the same geometric structure in the 



1 For some of the earlier attempts to quantizing the anisotropic spacetimes based on the old quantization methods 
see eg. [25M27I ] . See also Ref. [HI for a discussion on von-Neumann stability issues in anisotropic models in LQC. 

2 For some of the phenomenological implications see for eg. [43447| . 

3 In the loop quantization of perfect fluid there are several open issues such as that of deparametrization. See eg. 
Ref. [iil for a discussion. In this analysis, our treatment of perfect fluid will be at an effective Hamiltonian level. 
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pre-bounce and post-bounce phases? If not, then is such a transition across the bounce random, or 
does it follow a set of rules? If such rules exist, under what conditions are some transitions favored 
over others, and finally how do they depend on the matter content? 

In our presented analysis, based on effective spacetime description, we outline all the conditions 
under which particular geometric structures form in the Bianchi-I model in LQC. This result stands 
on an analog footing as the work of Jacobs' in the classical theory JH|. We show that there exists 
a hierarchy amongst different structures. Whether or not a certain geometric structure forms, 
depends on the ratio of shear and the expansion scalar. For large values of expansion normalized 
shear, cigar structures are more probable, whereas for small values of ratio of shear and expansion 
scalar, a point like or isotropic approach to the bounce occurs. Moreover, depending on the initial 
conditions, pre-bounce and post-bounce branches may have different geometric structures. That 
is, a universe which forms a cigar type structure in the contracting branch before the bounce, 
may form cigar, barrel, pancake or point like structure in the expanding branch (in the backward 
evolution). Interestingly, such a Kasner transition is not random. It follows selection rules which 
are determined by the initial relative strength of the matter and anisotropies, captured in terms of 
two anisotropy parameters. The Kasner transitions which take place across the bounce in Bianchi-I 
model in LQC remind us of the mixmaster behavior between Kasner phases in classical Bianchi-IX 
universe. However, there are important differences. As mentioned earlier, there are no Kasner 
transitions in the Bianchi-I model in GR. Thus, such transitions are purely quantum gravitational 
in origin. Further, unlike the Kasner transitions in the classical Bianchi-IX spacetime, there is only 
one transition across the bounce. And finally, the Kasner transitions in the classical Bianchi-IX 
model lie on the same branch (expanding or contracting) of the universe, whereas in LQC, they 
involve a non-singular transition between expanding and contracting branches. 

The plan of this work is as follows. In the following section, we start with the classical GR to 
provide a stage for discussion of the new results in LQC. We discuss the derivation of the primary 
dynamical equations in classical Bianchi-I spacetime using classical Hamiltonian constraint, and 
illustrate various geometric structures using the Kasner exponents. Here we provide a hierarchy 
among different geometric structures, which to our knowledge has not been discussed before in the 
classical theory. We discuss parameterization introduced by Jacobs [H] which proves very useful in 
the analysis of the formation of geometric structures. This is followed by a summary of the results 
in the case of stiff matter or the Zel'dovich solution, and dust and radiation. In Sec. Ill, we start 
with the effective Hamiltonian constraint of Bianchi-I model in LQC, and summarize the derivation 
of the effective dynamical equations which are then numerically solved. We analyze the cases of 
stiff matter, and dust and radiation separately, and classify all kinds of geometric structures and 
the allowed and forbidden transitions across the bounce. Defining a ratio in the parameter space, 
we also compute how favored a particular type of transition is for a given equation of state and 
the anisotropic shear scalar. We conclude with a summary of main results in Sec. IV. 

II. CLASSICAL THEORY: DYNAMICS, ANISOTROPIC HIERARCHY, JACOBS' 
PARAMETERIZATION AND EXAMPLES 

Since loop quantization is based on classical gravitational phase space variables, the Ashtekar- 
Barbero connection A\ and the triads Ef, in this section we start with a discussion on the classical 
theory in connection variables, and then relate it to conventional metric variables. In the first part of 
this section we obtain the dynamical equations resulting from the classical Hamiltonian constraint 
and discuss the Kasner exponents obtained from integration of the dynamical equations. We then 
use these Kasner exponents to elucidate the structure of the singularities in classical Bianchi-I 
model and establish a hierarchy amongst them using expansion normalized shear scalar. A useful 
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way to understand the properties of the anisotropic structure of the spatial geometry of spacetime 
is via parameters introduced by Jacobs jjj. After discussing these parameters, which we later use 
to classify the Kasner transitions, we illustrate some key properties of Bianchi-I spacetime with 
stiff matter (or massless scalar), and dust and radiation. 



A. Dynamical equations 



We consider a homogeneous and diagonal Bianchi-I spacetime with a manifold S x R where E 
refers to the flat spatial manifold having a R 3 topology. Due to the non-compact nature of the 
spatial manifold, one needs to introduce a fiducial cell V in order to define the symplectic structure. 
The fiducial cell is chosen to have a fiducial volume V = hhh where U are the coordinate lengths 
in the three spatial directions. The edges of the cell are chosen to lie along fiducial triads e l a . We 
denote the fiducial metric compatible with the co-triads uj l a as q a b- Using the underlying symmetries 
of the homogeneous spatial manifold, the Ashtekar-Barbero connection A l a and the triads Ef , can 
be symmetry reduced as: 



4 = e (k)- 1 &l and Ef = Pi k V- 1 , (2.1) 

where i = 1,2,3. The symmetry reduced connection c l and triad pj satisfy the following Poisson 
bracket relation 

{c\ Pj } = %vG 1 8) (2.2) 

where 7 = 0.2375 is the Barbero-Immirizi parameter fixed by the black-hole entropy calculation in 
LQG. The triads are related to the directional scale factors a* in the Bianchi-I metric 

ds 2 = -dt 2 + a\ dx 2 + a 2 2 dy 2 + o§ dz 2 (2.3) 

as 



Pi = hh\a2ds\, P2 = hk\aia 3 \, p 3 = hh\axa 2 \. (2.4) 
The classical Hamiltonian constraint for Bianchi-I spacetime in terms of the above connection 



and triad variables, for lapse iV = 1, is 27] 



^ci = 77 n T , {C1P1C2P2 + cyclic terms) + H ma , t t , (2.5) 
87rG7 z V 

where without any loss of generality we have chosen fiducial lengths to be unity, V = ^/piP2P3 
is the physical volume of the cell V and % m att denotes the matter part of the Hamiltonian. We 
consider the matter as a perfect fluid with a barotropic equation of state w = P/p, where P and 
p denote the pressure and the energy density of the matter content. 

The dynamical equations for the triad and the connection components can be derived using the 
Hamilton's equations: 

Pi = {puUd} = -8vrG 7 — A bi = {c h U cX } = 8irG 7 — ^ (2.6) 

oci dpi 

where a 'dot' refers to the derivative with respect to the proper time t. Using these equations, it 
is straight forward to obtain the classical evolution equations for Bianchi-I spacetime 

HiH 2 + H 2 H 3 + = 8nGp (2.7) 
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and 



H2 + H3 + Hi + Hl + Hl = SttGP (2.8) 



(and its cyclic permutations). Here Hi denote the directional Hubble rates, which are related to 
the time derivatives of the triad components, such as 

(2.9) 

a\ 2 \p 2 p 3 pij 

(and similarly for H2 and-f^). The energy density and pressure in the above equations, satisfy the 
conservation law 



p = -3H(p + P), (2.10) 

where H = (Hi + H2 + H%)/3 is the mean Hubble rate. Using the relation P = wp, the energy 
density can be expressed as follows 

p(t) = p a(t)- 3 ^ , (2.11) 

where p Q is the initial energy density, a = (aia2«3) 1//3 is the mean scale factor and a(t = 0) = 1. 

We now consider the parts of the expansion matrix in the Bianchi-I spacetime. For the diagonal 
Bianchi-I metric, the expansion tensor 9 a p is given by the covariant derivative of the unit fluid 
velocity, with its trace defining the expansion scalar 6 of a co- moving observer: 

6 = H 1 +H 2 + H 3 . (2.12) 

The traceless part of the expansion matrix defines the anisotropic shear which takes the following 
form for the diagonal Bianchi-I spacetime: 

Oij = diag (<7i, cr 2 , 03) (2.13) 

where Oi = Hi — ^6 are the anisotropy parameters. Since the shear matrix cr^- is traceless, by 
definition its components satisfy the following constraint 

ai + cr 2 + 03 = 0. (2.14) 

An important measure of the anisotropy of the spacetime is the shear scalar, defined as a 2 = 
a l:, (Tij, which can be expressed in terms of Hi as 

° 2 = E °i = I ((^l " H ^ 2 + - H 3 f + (H 3 - H,f) . (2.15) 

i 

In the classical theory, using the dynamical equations it follows that a 2 oc a -6 . Further, it 
satisfies the generalized Friedmann equation: 

H* = ^ P + \a 2 . (2.16) 
6 b 

In the absence of the matter, solution of the dynamical equations (|2.7p and (|2.8|) . gives the 
Kasner solution according to which the directional scale factors can be written as 

en oc t h (2.17) 
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where constants ki are the Kasner exponents. For vacuum Bianchi-I, these exponents satisfy the 
following two constraints 

ki + k 2 + k 3 = 1 and k\ + k\ + k\ = 1. (2.18) 



For general matter, the directional scale factors can not be expressed in terms of eq. (|2.17p with 
(12.1811 holding true. However, for the case of w = 1, it is possible to write similar constraints, given 

bTB 

kx + k 2 + k 3 = 1 and k\ + k\ + fc| = 1 - k 2 , (2.19) 

where k 2 is a constant which can be expressed in terms of the energy density and shear scalar 
as k 2 = | ( g^Jp/^f^/e ) • ^ or a mass l ess scalar field (for which the equation of state is w = 1) 
the constant k 2 can also be expressed in terms of its field momentum (Ps, which is a constant of 
motion) as- 



_ 2 / 4ttGP?/3 \ 
where, S 2 = \& 2 cP is also a constant motion. 

We now have the expression for k 2 in terms of the two constants of motion of the system for a 
massless scalar field. It is evident from the above equations that its value depends on the initial 
data and k 2 < 2/3. The Kasner exponents in terms of the directional Hubble rates can be given 
as ki = Hi/(3H). Utilizing this relation and Oi = Hi — H, it is straightforward to show that 
the constraints on the Kasner exponents given in eq. (I2.19P follow from eqs. (|2.14l and I2.15p . The 
second constraint in eq. ()2. 19j) can also be derived from the classical Hamiltonian constraint. 

The case of massless scalar field is interesting as it is the only matter, satisfying dominant 
energy condition, for which Bianchi-I universe does not isotropize as the universe expands. This 
is so because the energy density for w = 1 scales as p oc a~ 6 , in the same way as a 2 . Behavior of 
the energy density in relation to a 2 also shows that unless the matter is massless scalar field, the 
anisotropic shear always dominates p as the classical singularity is approached (for w < 1). Thus, 
for w < 1, a Bianchi-I spacetime near the singularity approaches the Kasner vacuum solution. 



B. Structure of the singularities and their hierarchy 

Presence of the anisotropic shear makes the structure of the spacetime near the singularity 
much richer in comparison to the isotropic spacetime. Recall that in an isotropic universe, the big 
bang singularity is always point like. However, in the Bianchi-I spacetime the geometrical nature 
of singularity depends on whether or not the scale factors in the three directions approach zero 
together, and can be classified as follows [1] 

• Point: The approach to singularity is said to be point like when all the scale factors tend 
towards zero simultaneously i.e. when all three Kasner exponents k\,k 2 ,k 3 > 0. These 
singularities occur when the energy density of the matter content plays an important role, 
and thus only occur for w > 1. At these singularities, the anisotropic shear does not play 
a dominant role. In the above sense, a point singularity is often also termed as isotropic 
singularity in Bianchi-I model. 
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Barrel: Barrel singularity occurs when one of the scale factors tends to a constant while 
the other two approach zero i.e. a\ — > const and 02,03 — > 0. The corresponding Kasner 
exponents for barrel type singularity are k\ = and k2,k% > 0. 

Pancake: Pancake singularity occurs when one of the scale factors tends to zero while the 
other two approach to constant i.e. a\ — > and 02, 03 —> const which means that the Kasner 
exponents k\ > and k2,k^ = 0. This singularity appears as one of the axisymmetric 
solutions for < w < 1 and does not take place for w = 1. 

Cigar: The approach to singularity is classified as cigar when one of the scale factors 
diverges while the other two approach to zero i.e. a\ — > 00 and 02,03 — > as the singularity 
is approached. For cigar type singularity one of the Kasner exponents is negative while the 
remaining two are positive i.e. k\ < and /C2, &3 > 



o 2 

1.0 r 



Cigar 

Pancake (=2/3) 



Point, Barrel, Cigar 



Barrel (=1/6) 



FIG. 1: Hierarchy of singularity structure in the classical theory based on a 2 /6 2 for Bianchi-I spacctimc. 



It is to be noted that the approach to singularity for the vacuum Bianchi-I spacetime is always 
either cigar or pancake like. Introducing the stiff matter (perfect fluid with w = 1) to the Bianchi-I 
spacetime brings in other types of the approach to singularity which were otherwise forbidden in 
the vacuum Kasner solution. The geometric structure of singularity is decided by the anisotropy 
present in the system. As remarked above, for a point like singularity to exist, the matter part 
should dominate over the anisotropy. On the other hand if the anisotropy wins over the matter, 
the point like approach to singularity becomes less likely. 

It turns out that one can define a hierarchy of the various approaches to the singularities based 
on the value of a 2 /6 2 , the normalized dimensionless shear scalar. Using the definitions of the shear 
and the expansion scalar, the normalized shear scalar can be expressed in terms of the Kasner 
exponents k{ as follows: 

a 2 2 

= -{l-3{k 1 k 2 + k 2 k 3 + k 3 k 1 )). (2.21) 

The ratio a 2 /9 2 takes the smallest value when all three Kasner exponents are positive which is the 
case for point singularity. It is also evident from the expression for a 2 /9 2 and the definitions of 
the various singularities that for a pancake singularity a 2 /9 2 = 2/3, for barrel singularity 1/6 < 
a 2 /9 2 < 2/3 and for cigar singularity a 2 /9 2 > 1/6. This hierarchy of the singularities based on the 
normalized shear scalar is shown in Fig.[TJ From here we can say that the point like singularity is the 
least anisotropic structure of singularity while cigar corresponds to highly anisotropic situations. 
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C. Jacobs' parameterization 



It is convenient to capture the information contained in these, by using parameters introduced 
by Jacobs Ref. One first defines two parameters e and ip as: 



4irGp 



(2.22) 



sin (V0 = \/j?4, sin^ + ^=V|4 ^(^ + ^-)=\H- (2-23) 



and 



2cr 2 ' V 3 7 V 2cr 2 ' V 3 / V 2cr 

We will refer to these transformed parameters as Jacobs' anisotropy parameters in the rest of 
the article. The ranges of the Jacobs' anisotropy parameters are |e| = (0, oo) and ip = [0, 2tt/3\. 
Note that in the classical theory e = is reached only in the isotropic limit of the model. Similarly, 
e = oo is also reached only for a vacuum Bianchi-I model or when the approach to singularity is 
dominated by anisotropics. 4 

It is evident from the eqs. (|2.22l and l2.23p that e refers to the ratio of the total anisotropic shear 
present in the spacetime to the energy density while the parameter ip tells how the total anisotropy 
is distributed in the three directions. Thus, one has complete information about the shear matrix 
Oij of the spacetime, when both e and ip are specified. 

The classical dynamical equations in GR can be analytically integrated for constant equations of 
state w. For < w < 1, upon integration, the Einstein's equations give the following expressions 
for the scale factors in terms of Jacobs' parameters [H] 

2 sgnfe)^ 



, 4a SM+^'V 3(l-») 

alt) ■} — — 1/2 (0 < w < 1) 

1 4a(f) 3 C 1 -'»>+e a 1/8 -|e| ' 



(2.24) 



a(t) l-2 ^(4+^-1/2 {w = 1) 



where 



Zi = sin (ip) , Z 2 = sin fip + , Z 3 = sin + . 
Also, the expressions for the directional Hubble rates are given via, 

= [ Wp^I {0 - w - 1] (2 - 25) 

where p Q is the initial energy density. In the following subsections we discuss the properties of 
these solutions for the case of massless scalar, dust and radiation. 



4 In the stiff matter universe, a non-zero finite |e| in the initial data ensures that |e| will never be zero or infinity 
during the evolution, not even if a 2 and p diverge as the singularity is approached. In LQC, since the energy 
density and shear both are finite during the entire evolution for Bianchi-I spacetime, e = is possible when a 2 = 0, 
which is the isotropic limit of the model and e = oo takes place only when p — which is only allowed for the 
vacuum case. 
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D. Stiff matter 



The solution of Bianchi-I spacetime with stiff matter is also referred to as Zel'dovich solution [5C 
For such a universe where w = 1 the energy density varies with the mean scale factor as p oc a~ 6 , i.e. 
in the same way as the shear scalar a 2 in the classical theory. Thus, the ratio of the shear scalar to 
the energy density is a constant, which means that the Bianchi-I spacetime filled with stiff matter 
never isotropizes. Further, depending on the initial conditions, it is also possible that the matter 
density dominates over anisotropy and the approach to singularity is point like. Hence, such a 
universe does not, in general, behave like vacuum Kasner universe as the singularity is approached. 
For a Bianchi-I universe with stiff matter, |e| is a constant of motion in the classical theory, and 
the Kasner like behavior occurs only when the anisotropic contribution to the spacetime curvature 
dominates over the contribution due to the matter. This requires e to be greater than a critical 
value e*, where for any e < e*, the approach to singularity will be always point like and the Kasner 
behavior will be absent. We now find out this critical value e. 

We first make the following transformation of the parameter e 0]: 

(2.26) 



This transformation restricts the range of new anisotropy parameter 8 to < \5\ < 1. In terms of 
5 and ip, the expressions for the directional scale factors and Hubble rates for stiff matter universe 
take a much simpler form, 

a l (t) = a lo a(t) 1 ±^ J and H&) = ^t5^ • ^ 

Here aj is the directional scale factor at the initial time. From the above equations, it is clear 
that for all of the Hubble rates to be positive, which is the requirement for a point like approach 
to singularity (in the expanding branch), 5 must be given by 1 ± 2\5\ > 0. Thus, for \6\ < 1/2 the 
approach to singularity will always be point like, irrespective of the value of ip. As we will discuss 
in the next section, for 5 > 1/2, the approach to the classical singularity can be point like, barrel 
or cigar. Using the definition of 5, we find that the critical value of e below which the approach to 
the classical singularity is always point like turns out to be e* = 2/\/3. Existence of such a critical 
value is an important property of the Bianchi-I spacetime with stiff matter, and we will see in the 
next subsection that such a critical value of anisotropy parameter is absent for uu < 1 for which 
the point like singularity is forbidden. Integrating Einstein's equations, the mean scale factor for 
w = 1 universe can be written as a function of time as 

a(t) oc (t + t s ) l/3 (2.28) 

where t s is an integration constant. Now let us write the classical directional scale factors in terms 
of the Kasner exponents ki as follows: 

Oi(t) oc (t + t s ) k > (2.29) 

where ki are the Kasner exponents. The anisotropy parameters 5 and ip can now be expressed in 
terms of the directional scale factors and the Kasner exponents as follows - 



9 



and 



Am 3(fci - fc 2 ) 

<>W = / ■ ■ ,ry » • (2.31) 

2 (sm(Zi) - sm(Z 2 )) 

Bianchi-I spacetime with stiff matter admits all kind of singularity except the pancake type, as 
the Kasner exponents for the pancake singularity are k\ = 1 , k 2 , ks = and they do not satisfy 
the constraint on the Kasner exponent for w = 1 given in eqs. (j2. 19|) . The axisymmetric solution 
could be either cigar or barrel depending on the value of 5 and ip. We notice that if the anisotropy 
parameter (8) is small enough, the only possible singularity is point. For higher values of 6, 
barrel and cigar singularity form depending on the value of the parameter ip. The various possible 
structure of singularity in the classical Bianchi-I model for stiff matter universe are shown in the 
second column of Tables U and [TT] in the following section. 

E. Dust (w = 0) and radiation (w — 1/3) 

Solution for a general equation of state in the range < w < 1 is given via eqs. (|2.24l and l2.25|) . 
The energy density depends on the mean scale factor as p oc a(t)~ 3 ( 1+w ^ and hence the ratio of the 
anisotropic shear scalar to the energy density is given via 



2 

17 oc a(t)-^ w ^ . (2.32) 
P 



It is evident from the above equation that the shear always dominates over the energy density 
for w < 1 as the singularity is approached i.e. a(t) — > which implies that the point like singularity 
does not occur. Near the singularity the spacetime behaves similarly to vacuum Bianchi-I given 
via the vacuum Kasner solution. Moreover, the axisymmetric case is either a cigar or a pancake 
singularity, and the barrel singularity does not form. 5 In the asymptotic limit when a — > oo, it is the 
energy density which is more prominent and contribution of the anisotropic shear to the extrinsic 
curvature become negligible compared to that of the energy density. Hence the isotropization takes 
place. This is an important distinction from w = 1 where the universe never isotropizes, not even 
in the asymptotic limit. It can also be concluded from the eq. (|2.32p that smaller the equation 
of state, faster will be the isotropization. We now discuss some of the main features of the dust 
filled Bianchi-I universe. (The features for the radiation filled case can be deduced similarly). 
The behavior of the directional scale factors for a dust filled universe (w = 0), also called the 
Heckmann- Shucking solution, is given as [5l| 

Oi(t) = a l0 t k *(t + t f 3 - k > (2.33) 

where t D is the initial time, and k{ are the vacuum Kasner exponents satisfying the eqs. (|2.18p . As 
the classical singularity is approached i.e. when t then it is straightforward to see from the 

eq. (12.33h that the directional scale factors behave as a,i(t) oc t ki .® This precisely is the Kasner 



5 For barrel singularity, ki = 0, k-z, &3 > 0. These values of the Kasner exponents do not satisfy the constraint for 
the vacuum Bianchi-I given by eqs. (|2.18[) therefore barrel singularity does not form for dust and radiation filled 



Bianchi-I spacetime. 
For t < U 
expression. 



6 For t to, t + to ~ t which is a constant and therefore the time dependence comes from the t ki part of the 
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vacuum solution to the Bianchi-I spacetime. This shows that, near singularity, w = universe 
behaves like vacuum solution to Bianchi-I. Note that the Kasner solution is only a limit which is 
approached when t <C t Q , the universe never actually becomes an exact Kasner vacuum universe. 
On the other hand the asymptotic limit t S> t gives the following expression for the directional 
scale factors ai(t)dust|t^oo oc i 3//2 and in this limit the isotropization takes place. Following similar 
argument for a radiation filled universe, one finds that the spacetime behaves like vacuum Kasner 
universe as the singularity is approached while in the asymptotic limit the spacetime isotropizes and 
the directional scale factors behave as aj(i) ra d|t->oo oc i 1 / 2 . It is worth noting that the asymptotic 
behavior of these scale factors for dust and radiation filled Bianchi-I spacetime is the same as that 
of the scale factor of an isotropic spacetime filled with dust and radiation respectively. 



III. EFFECTIVE DYNAMICS IN LQC AND KASNER TRANSITIONS 

In the previous section, we discussed the way various geometric structures arise in the Bianchi-I 
spacetime with stiff matter, dust and radiation in the classical theory. In this section, we will 
examine the formation of these structures in LQC using effective dynamics. We will study the 
structure of the spatial geometry on both sides of the bounce by computing the Kasner exponents 
and investigate the various possible transitions as the spacetime transits from an expanding to a 
contracting branch during the backward evolution of given initial data. The dependence of these 
transitions on the matter content and the anisotropy will also be explored in detail. We will, first, 
discuss the effective dynamics of Bianchi-I spacetime, and then in the following subsections we 
will present the results of numerical solution of the dynamical equations for stiff matter, dust and 
radiation filled universes. 



A. Effective dynamics 

The loop quantization of Bianchi-I is performed by expressing the classical Hamiltonian in 
terms of the elementary variables of LQG, namely the holonomies of the connection A l a taken 
over the closed loops and the fluxes of the triads Ef and promoting them to their correspond- 
ing quantum operators. This results in discrete quantum difference equations which govern the 
evolution of the spacetime. The discrete difference equations hence obtained are manifestation of 
the discrete quantum geometry inherent in loop quantization. LQC also provides an approximate 
effective description of the underlying quantum geometry of the spacetime, known as the effective 
Hamiltonian description. The effective description is obtained through a faithful embedding of the 
finite dimensional classical phase space in the infinite dimensional quantum phase space of LQC 
based on the geometric formulation of quantum mechanics (5^, [H^]. Through a proper choice of 
semiciassicai states, this procedure furnishes an effective Hamiltonian which contains the leading 
order corrections due to the discrete quantum geometry. Using the Hamiltonian hence obtained, 
the effective equations of motion can be derived. These modified equations of motion introduce 
quantum geometric corrections to the classical dynamics of the spacetime and yield the modified 
dynamical equations. The dynamics hence obtained has been shown to be in excellent agreement 
with the full quantum evolution of semiciassicai states using numerical simulations in many situa- 
tions including in the presence of anisotropies lol-l^. El, EH El, [II |33( |. Our analysis will assume 



validity of the effective Hamiltonian approach for the Bianchi-I model. 

The effective Hamiltonian constraint for Bianchi-I with choice of lapse N = 1 can be written as 
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H, S3, S3 7 



1, 1/ ,17 1 /sin (/^ci) sin (// 2 C2) \ , ^ 

«eff = rtgrav + «matt = ~ p „ 9T , = = P1P2 + cyclic terms + 7i matt (3.1) 

where % ma tt is the matter Hamiltonian and /Zj are given via 



/Ui = A,/^-; fi2 = K—^- and fi 3 = X A — — . (3.2) 
V P2P3 V PiP3 V 

Here A 2 = 4^7^ arise from the discrete quantum geometry of the underlying spacetime pre- 
dicted by LQG. The dependence of /2j on the pi is unique in the sense that any other choice of 
the function fii(pi) would lead the physical results to depend on the rescaling of the edges of the 



fiducial cell and the shape of the fiducial cell would affect the physical results [44j]. 

The equations of motion of the triads and connection can be computed using the Hamilton's 
equations of motion- 



Pi _ _ _ 

Pi = — r (sin (/2 2 c 2 ) + sin(/2 3 c 3 )) cos (fiiCi) (3.3) 
7A 



and 

1 

ci = 



C2P2 cos (// 2 c 2 ) (sin(/uici) + sin(/x 3 c 3 )) + c 3 p 3 cos (p 3 c 3 ) (sin(/iici) + sin(/y, 2 c 2 )) 

P17A [ 

- cipi cos (fiici) (sin(/2 2 c 2 ) + sin(/2 3 c 3 )) - 2 ^ 2P3 f s in(/2 2 c 2 ) sin(/i 3 c 3 ) 

+ sin(/i 3 c 3 ) sin(/2ici) + sin(/2ici) sin(/2 2 c 2 ) J + 8irG ^^t ma,tt . (3.4) 

/J "Pi 

Similarly one can obtain equations for p 2 , p 3 , c 2 andc 3 , and also Sj/aj [jj ]. 

The values of the directional triads and the connection at any time can be computed through 
numerical evolution of initial data specified at a given instant of time. The directional Hubble 
rates, the expansion scalar and the shear scalar can be computed using the definitions given in 
eqs. (|2.9l 12.121 and I2.15P respectively. The dynamical equations of motion for pjandq (given in 
eqs. 13.31 and 13.41 and their cyclic permutation) modify the Einstein's field equations. From the eq. 
(|3.3p it is evident that since the maximum value of sin(0) is 1, the value oipi/pi never diverges and 
hence the modified effective equations provide a singularity free solutions, causing bounces of the 
directional scale factors (or the bounce of the mean scale factor). The energy density, expansion 



scalar and the shear scalar are also bounded by their respective upper maxima 35(. Further, all 
strong singularities for the matter with equation of state w > — 1 have been shown to be resolved 
0- 

Before we go in to detailed analysis for different matter components, it is important to make 
following remarks. 



7 The effective Hamiltonian considered here does not include 'inverse volume corrections' which are sometimes 
included in studies of spatially compact models in LQC. In our analysis the Bianchi-I model is considered with 
a non-compact topology. Even if one takes a compact model, since the physical volume the fiducial cells in the 
simulations presented here are larger than Planck volume, the corrections introduced due to inverse triad terms 
would remain negligible. 
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Remark 1: In the classical theory, depending on the matter content and the anisotropy, Bianchi-I 
spacetime demonstrates various geometrical structures of singularity; for example if one of the 
three directional scale factors diverges and the other two approach to zero as the singularity is 
approached, then the structure of the spatial geometry of the spacetime tends to an 'infinite' 
cigar. In LQC, for the same initial data, the structure near the bounce is rather a 'finite' cigar. 
This is so because in LQC, due to the nonsingular bounce the singularity is avoided and all the 
scale factors take non-zero and finite value near the bounce. Similarly, if the classical trajectory 
for a given initial data collapses to a barrel, point or pancake singularity, then in LQC, these 
geometric structures form for the same initial data, but are finite. Although LQC and GR 
trajectories are in agreement with each other in the low curvature regime and an observer sees 
both trajectories tending towards the same type of geometrical structure as the high curvature 
regime is approached, the actual spatial geometry near the bounce in LQC is non-vanishing and 
is a finite version of the geometry in GR near singularity. 

Remark 2: In the following analysis we always specify the initial data in the expanding branch 
(post-bounce) of the mean scale factor and perform a backward time evolution through the bounce 
to the contracting branch (pre-bounce). A Kasner transition in the following is defined as the 
change in geometric structure from the post-bounce phase of the mean scale factor a, as observed 
by an observer approaching the bounce, to the one in the pre-bounce phase as will be observed by 
the same observer if she/he evolves forward in time in the contracting branch. Thus, even though 
the evolution occurs in only one direction (i.e. backward), the geometric structure in the pre- and 
post-bounce phases is defined in the approach to the bounce. Note that the nature of transition 
is unaffected by providing the initial data in the contracting branch of the mean scale factor and 
evolving in future time direction. 

B. Bianchi-I LQC with stiff matter 

We perform several numerical simulations (over 100 simulations for each range of 6) for various 
anisotropy parameters to solve the system of differential equations (13.31 and I3.4p for a variety of 
initial data. We use 4th order adaptive Runge-Kutta for the numerical solutions to the equations. 
Then, for various ranges of the anisotropy parameters we find out the structure of spacetime by 
calculating the values of the Kasner exponents. This way a wide class of transition of structures is 
obtained, based on the initial anisotropy. It also turns out that the structures of spacetime on two 
sides of the bounce are not the same although the anisotropic shear parameter S 2 is conserved on 
the two sides of the bounce in the asymptotic limit. This provides a rich and interesting transitions 
of structures across the bounce. Fig. [2] shows an example of a Kasner transition across the non- 
singular bounce. We find that in the post-bounce phase, two scale factors are decreasing and one 
is increasing as the bounce of the mean scale factor is approached, i.e. a cigar is formed. Similarly, 
in the pre-bounce phase, two scale factors decrease while one increases, on approach to the bounce 
in the forward time evolution. Thus, again a finite cigar forms. Hence, this figure corresponds to 
a cigar-cigar transition. (In all the plots, presented in the following, structures are defined in this 
same way). 

In all numerical simulations presented here, the initial data is given at t = without any loss 
of generality and the classical singularity occurs at t = —t s . The value of t s depends on the initial 
data. The initial data set consists of the values of the connections (ci, C2), the triads (pi,P2iP3), the 
energy density p at time t = 0. The value of the remaining connection component C3 is computed 
by satisfying the effective Hamiltonian constraint given by eq. (|3.ip . The parameters tp and 5 are 
computed via eq. (12.301 and I2.31|) by using the Kasner exponents and the directional scale factors, 
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t 



FIG. 2: The left plot shows the cigar-cigar transition of structure of the spatial geometry of the spacetime 
for l/v3 < S < 1, the dotted lines represent the corresponding classical trajectory and the right plot shows 
the evolution of the mean scale factor. The spacetime undergoes a cigar-cigar transition while the mean 
scale factor bounces and remains non-singular throughout the evolution. 




FIG. 3: This plot shows the Hamiltonian constraint during the evolution of the effective equations. 
I'Hcff/'Hmattl and I'Hcff/'Hgravl remain sufficiently small at all the time which indicates the high numeri- 
cal accuracy of the results. 

which can be easily obtained from the triads pj. This provides the complete specification of the 
initial condition of the system. Also, the constant ip c (which will be used further in this section) is 
given in terms of the anisotropy parameter (<5) as ip c = sin -1 (|<5|/2). 

Fig. [3] shows an important feature of the numerical simulations, the absolute value of the ratios 
of effective Hamiltonian to the matter Hamiltonian and gravitational part of the Hamiltonian are 
of the order of 10 -8 or less during the evolution, which implies that the effective Hamiltonian 
constraint is satisfied with such a high order of accuracy. This, in turn, assures the numerical 
accuracy of the results presented in this article up to the order of 10~ 8 . We consider such a 
numerical accuracy to be good enough for our analysis because the physical quantities of interest 
are of order much higher than 10 -8 and hence the numerical errors can safely be neglected. 

We now discuss the cases of various values of the anisotropy parameter 5. 

1. 1/V3 < 5 < 1 

This range of anisotropy parameter 5 corresponds to high anisotropy and it is the only range in 
positive 5 region when cigar to cigar transition is allowed. All possible transitions depending on tjj 
are shown in the Table-I. Out of the total range of ip = [0, 2tt/3] , cigar to cigar, of all the transitions, 
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occupies the biggest range. As we will see in the following subsections that cigar-cigar transition is 
forbidden for \5\ < l/v3. This signals that cigar-cigar transition is a highly anisotropic transition. 
Axisymmetric solutions exist for ip = — ip c and ^ + ip c ) on the expanding side of the universe 
and for Tp = (| — ip c a,ndip c ) on the contracting side. As discussed in Sec. II, the parameter 5 is 
related to the ratio of the shear scalar to the energy density of the matter field and larger value of 
5 means that the anisotropy dominates over the matter energy density. This is partially the reason 
why there are no point like to point like transitions allowed for this range of 5. 



{(A=- -^,5 = 0.8} 
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Barrel 






Cigar 
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{^ = ^-^,5 = 0.8) 
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FIG. 4: These plots show various transitions of structure across the bounce depending on various values of 
ip for I/a/3 < 5 < 1. 



{if/ = 0.58643 1. 6 = 0.8. if/ c = 0.675 132) 



Point 



Cigar 



FIG. 5: This figure shows the transition of Kasner exponents (from cigar to point like) zoomed in near the 
bounce for S = 0.8 and ip = 0.58543. 



Various possible transitions for a choice of 5 = 0.8 are shown in Fig. [31 The transition for 
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5 = 0.8 and tp = 0.586431 (which correspond to cigar to point like transition) zoomed in near 
the bounce has been shown in Fig. [5] which clearly portrays the rich dynamics near the bounce 
and resulting transition in the Kasner exponents. It is straightforward to see that the values 
of the Kasner components change across the bounce leading to the corresponding transition in 
the structure of the spatial geometry. Similar to the Kasner transitions in classical mixmaster 
phenomena, one of the Kasner components retains its behavior and the remaining two switch their 
role or two of them retain their behavior while the remaining one changes its sign. In other words, 
if k\ < 0and/c2,fc3 > (cigar) in the expanding branch, then after the transition takes place the 
Kasner components are either k\,k2 > OandA/j < (as shown in the cigar-cigar transition where 
one of the Kasner exponents remains positive on both sides while the remaining two change their 
sign) or k\ > 0and/c2,&3 > (as shown for the cigar to point like transition where only of the 
exponents changes the sign while the other tow remain positive across the bounce). The table also 
shows that if there is an axisymmetric barrel structure on one side, it goes to cigar type structure 
across the bounce. There are no barrel to barrel, or barrel to point like transitions allowed for 
this range of 5. The transitions involving barrel happen for a certain value of tp, whereas other 
transitions have a range of tp associated to them. 

2. H=l/>/3 

5=1/ \/3 is very interesting case in the sense that this is the only case where the barrel 
to barrel transition occurs. In all the other cases if there is barrel structure on one side of the 
bounce, it becomes either cigar or point like on the other side. Except for the two values of ip 
when the barrel structure is allowed, for all the other ip the transitions are the point like to cigar 
transitions. Moreover, the cigar-cigar transition is forbidden in this range. The barrel to barrel 
transition occurs at particular values of ip unlike for other transitions for which there are ranges 
of ip. Fig. [6] shows the Kasner transition for 5 = 1/V~3. Fig. [7] shows one of the transitions of 
Fig.[6] {5 = 1/^/3, V> = 2tt/3 - 0.1) zoomed in near the bounce. In a barrel-barrel transition one 
of the Kasner exponents remain the same while the other two interchange their value across the 
bounce. It is important to note that the transitions involving barrel occur at a fixed value of tp for 
a given 5 while other structures happen for a range of values of ip. This feature is shown in Fig. [6] 
- the barrel-barrel transition occurs at tp = 2-7r/3, and if ip is decreased by a small amount, the 
structures on both sides of the bounce change and the resulting transition becomes point like to 
cigar. Similarly, when tp increased slightly, cigar to point like transition turns up. 

3. l/2<(5<l/ v / 3 

In this range of 5, point like to point like transition appears but there are no cigar-cigar tran- 
sitions anymore. Also the barrel-barrel transition are now forbidden. The absence of cigar-cigar 
transitions and the appearance of point like to point like transition are consistent with the fact that 
the anisotropy is now smaller than the previous cases, and the matter dominates over the anisotropic 
shear. Fig. [8] shows various transitions for 5 = 0.53. One of the transitions (5 = 0.53, ip = 0.04188), 
which correspond to point like to point like transition, has been zoomed in near the bounce in 
Fig. EE It is worth noting that barrel to point like transition occurs at a fixed value of -0 = ip c — n/3 
and if ip is changed even for a small amount, the transition of structures changes. As shown in 
figure if ip is decreased, the transition is point like to cigar and when increased it becomes cigar to 
point like. Again, at ip = 2ir/3 — ip c the transition involves the barrel structure in the expanding 
branch. 
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FIG. 6: Various transitions of structure across the bounce for 3=1/ V3. The Barrel-Barrel transition occurs 
at ip = 2ir/3, a small deviation from ip = 2it/3 results in other transitions for example when -0 decreased, 
the transition is Point-Cigar and when increased, the transition is Cigar-Point. 
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FIG. 7: This figure shows the transition of Kasner exponents (from point like to cigar) zoomed in near the 
bounce for S = 1/V3 and if) = 2ir/3 — 0.1. 



4- 1*1 < 1/2 

This range of 5 is highly isotropic range in the sense that the matter dominates over the 
anisotropy so much that the transition is always point like to point like except for the two special 
cases namely (5 = 1/2, ip = 7r/2)and(<5 = —1/2, ip = 7r/6), for which the transition is point like 
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Table-I 





Expanding branch 


Contracting branch 


l/v / 3< S < 1 


Cigar 

o < v < if - 


Cigar < V < f - i>c 

Barrel ip = | — ip c 
Point | - ip c < ip < ip c 

Barrel ip = ip c 
Cigar Vc < ip < x - V'c 


Barrel 

V> = f ~ ^ 


Cigar 

V> = IT - 


Point 


Cigar 

^ - V c < V> < f + i> c 


Barrel 

= § + 


Cigar 

^ = f + A 


Cigar 

f + xj: c < V < ^ 


Cigar 

f + iP c < ip < f 


<y = i/V§ 


Barrel ^ = 


Barrel -0 = 


Cigar < ip < § 


Point < ip < § 


Barrel V = f 


Barrel ip = ^ 


Point f < V < 


Cigar f < V < ir 


1/2 < (5 < 1/V3 


Point 

o < v < V-c - f 


Point 

o < V < ^ - | 


Barrel 

$ = A - f 


Point 

V» = V'c - f 


Cigar 

v> c - f < v> < ¥ - ^ 


Point 


Barrel 


Point 

(iP = f _ ^) 


Point 


Point (2jC - C < i> < ip c ) 

Barrel (ip — ip c ) 
Cigar (ipc < tp < 7T - ip c ) 

Barrel (ip = tt — tp c ) 
Point (tt - < ^ < if) 



TABLE I: Transition of structures across the bounce for (1/2 < 6 < 1) 



to barrel. The barrel to point like transitions again appear at a fixed value of ip. This feature is 
shown in Fig.llOi In Fig. [TT] the transitions for 5 = 1/2, and^ = 7r/4 has been zoomed in near the 
bounce. The transition involving barrel is the axisymmetric case and all the other transitions are 
always point like to point like transitions. For \5\ < 1/2 the only structure possible is a point like 
and therefore the transition for \5\ < 1/2 is always point like to point like. As one decreases \5\ 
further towards zero, the universe becomes more and more isotropic. 

Tables HI and [TT1 summarize the geometric structures on both sides of the bounce for various values 
of 5 in the range —1 < 5 < 1. One can see that as the value of |<5| increases, more anisotropic 
transitions start showing up i.e. for \5\ < 1/2 only point like to point like transitions are allowed. 
For \S\ > l/\/3 most prominent transition is cigar-cigar. In this way, there is a preference of 
certain transitions over the others depending on the value of |<5|. Table- HTTI shows all the allowed 
and forbidden transitions for various values of the anisotropy parameter \6\. Since, there are 
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FIG. 8: This figure shows the transition of structure across the bounce for 1/2 < 5 < when S — 0.53. 

= 0.0418879, 8 = 0.53, ip c = 1 .23273) 



Point 



Point 



FIG. 9: This figure shows the transition of Kasner exponents (from point like to point like) zoomed in near 
the bounce for 5 = 0.53 and ip = 0.04188. 

favored sets of transitions for different values of 5, let us now define the "transition fraction" T of 
a transition as the ratio of the range of ip for which a particular transition takes place to the total 
range of ifj- F° r example, the point like to point like transition takes place for 1/2 < |<5| < l/\/3 
for the following values of ip: 



2?r , 

+ ( — - (7T - Vc) 



IpP^P = 

Therefore the transition fraction of point like to point like transition is 



(3.5) 
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FIG. 10: The left figure shows the axisymmetric transition for S = 0.5 and right figure shows the transition 
for the same value of S and ip = ?r/4 which gives the point like to point like transition. 




FIG. 11: This figure shows the transition of Kasner exponents (from point like to point like) zoomed in near 
the bounce for S = 1/V3 and ip — 2n/3 — 0.1. 



T P 



4 (ip c - tt/3) 
2vr/3 



For \5\ < 1/2 the transition fraction of point like to point like transition is 1. Similarly, T for 
various transitions can be calculated: 

(i) for point like to point like transition 



- | sin 

1 





- 1 ' -L 



(1/2 < \6\ < 1/V3) 

m < i/2) 

(otherwise) 



(3.6) 



(ii) for point like to cigar transition 



T P 



§U-2sin- 1 f^ 



(1/V3 < \5\ < 1) 
(1/2 < \S\ < I/VS) 
(otherwise) 



(3.7) 



(hi) and, for cigar-cigar transition 
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Table-II 





Expanding branch 


Contracting branch 


1*1 = 1/2 


Point 

1> ? I (5 > 0) 


Point V ^ | 
Barrel xj) = \ 




Point 

V> 7^ f (S < 0) 


Point V + § 
Barrel i/j = ^ 


0< \S\ < 1/2 


Point (V0) 


Point (V0) 


-1/V3 < <5 < -1/2 


Point 

o < V < V> c 


Point < -0 < - f 
Barrel ip = ipc — § 
Cigar - f < V> < if ~ ^ 
Barrel tp = — ip c 
Point 3f - -00 < "0 < "0c 




Barrel 

^ = ^c 


Point 

"0 = "0c 




Cigar 

0c < < 7T - -0 C 


Point 

•0c < if) < 7T - -0 C 




Barrel 

-0 = 7T - -0c 


Point 

1p = TT-1p c 




Point 

7T - -0c < "0 < -X 


Point 

7T - "0 C < "0 < — 


5 = -l/V3 


Barrel (-0 = 0) 


Barrel (-0 = 0) 




Point (0 < V < f ) 


Cigar (0 < ip < f ) 




Barrel (0 = § ) 


Barrel (0 = § ) 




Cigar (f < -0 < 3f ) 


Point (f < -0 < it) 


-1 < 5 < -l/\/3 


Cigar 

o < v < 5 - V'c 


Cigar 

o < v < f - V'c 




Barrel 

V> = f - Vc 


Cigar 

i/> = I - i>c 




Point 

f - -0 C < -0 < Vc 


Cigar 

f - -0c < V < 




Barrel 

"0 = "0c 


Cigar 

"0 = 4>c 




Cigar 

Ik < V> < x 


Cigar Vc < ip < x - 0c 
Barrel = ^ - -0 C 
Point ^-<0 c <<0<§+<0c 
Barrel -0 = ? + -0 C 
Cigar f + -0c < ^ < — 



TABLE II: Transition of structures across the bounce for (— 1 < 6 < 1/2) 



w = p(i-™ '(*)) (W3<|*|<i) (38) 

[ (otherwise) 

Fig. [12] shows the way transition fraction of point like to point like transition depends on \5\. 
T for point like to point like transitions is highest for smaller values of 5. As |<5| increases beyond 
1/2, the transition fraction of point like to point like transition decreases and cigar to point like 
transition starts increasing and attains a maximum. At this maximum, point like to point like 
transition has zero transition fraction and as \5\ is increased beyond this value, the cigar-cigar 
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FIG. 12: Fraction of various transitions. The (red) solid line corresponds to point like to point like transition, 
(blue) dotted line represents the cigar to point like transitions and the (orange) dashed line shows the 
likelihood of the cigar-cigar transition. The (blue) thick line for 0.5 < \6\ < 0.57 shows that the transitions 
involving barrel take place in this range, however the transition fraction T corresponds to a set of measure 
zero for such transitions. 
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TABLE III: This table shows that depending on 5, some transitions are favored over others. As \S\ increases, 
the transitions across the bounce involve more anisotropic structures. For \S\ < 1/2 only transition possible 
is point like to point like and as \S\ increases more anisotropic transitions take place and for \S\ > l/y/3 the 
point like to point like transition is forbidden. 



transitions begin to increase from zero while T for cigar to point like transitions decrease. As the 
anisotropy increases, cigar-cigar transitions become more prominent than the rest of the transitions. 

Remark 3: As noted in this section, barrel structures are special in the sense that they occur at 
a particular value of tp for a given 5 and so do the cigar-barrel or point like to barrel transitions. 
According to the definition of T, the transition fraction for any transition involving barrel is 
essentially negligible (corresponds to a set of measure zero). Table- Hill shows that transition to 
(from) barrel is possible only for 1/2 < |<5| < l/y/3. The range of transitions involving a barrel is 
marked in the Fig. [12] by a (blue) thick line. 

Remark 4: Tables U and [TT] show an important symmetry regarding the value of 5 and structure 
of the spatial geometry near the bounce. A positive (negative) 5 in the expanding branch forms 
the same structure, for a given value of ip, as negative (positive) 5 in the contracting branch. For 
example, cigar structure is formed for < tp < (2ir/3 — ip c ) in both the cases when 1/V~3 < 5 < 1 
in the expanding branch and — 1 < 5 < —l/y/3 in the contracting branch. This behavior can 
be understood by writing 5 in terms of the mean Hubble rate H as 5 = ±i \f o 2 /Q. As the 
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universe makes a transition from expanding to contracting branch, H changes from positive to 
negative. This leads 5 to change it sign across the bounce. For example, if one starts with a 
positive (negative) value of 5 in the expanding branch, then, during the backward evolution across 
the bounce, 5 becomes negative (positive) in the contracting branch. 



C. Dust (w = 0) and Radiation (w = 1/3) 



In this subsection we discuss the Kasner like transitions for the equation of state w = and w = 
1/3. From eq. (|2.32j) . it is clear that the behavior of anisotropy with respect to the matter energy 
density is different in the two limits i.e. a — > and a — > oo. Recall that in the classical theory, 
the shear scalar always dominates over the energy density (for w < 1) and the universe behaves 
like a vacuum Bianchi-I near the singularity, whereas in the asymptotic limit universe isotropizes 
giving the directional scale factors to be ai(t) oc i 3//2 for dust (w = 0) and aj(£) oc t 1 / 2 for radiation 
(w = 1/3). In LQC, as the expanding branch of the spacetime is evolved backwards in time the 
bounce of the mean scale factor happens well before the vacuum solutions are reached. Since 
there are no point like or barrel structures allowed for the dust and radiation filled universe, the 
transitions involve only pancake and cigar structures. Also, since the classical behavior of dust and 
radiation filled universe near the singularity is the same (i.e. the vacuum Kasner solution), they 
allow similar types of transitions across the bounce. Table IIVI shows all possible and forbidden 
transitions for dust and radiation filled Bianchi-I universe in LQC. It also turns out that other 
equations of state in the range < w < 1 follow the same selection rule given in the table. We 
find that in all the transitions the structure is cigar on at least one side of the bounce. There is 
no pancake-pancake transition present. 
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TABLE IV: Transitions for < w < 1. (P— sToint, B-KBarrel, Pc-^Pancake and C— s-Cigar) 



For dust and radiation models, the pancake type structure comes out to be special (similar to 
barrel in the stiff matter universe) in the sense that it appears only as an axisymmetric solution 
occurring only for certain values of ip (ip = 7r/6 when e > and ip = tt/2 when e < [5!]) for a given 
|e|. Moreover, pancake to pancake like transition is forbidden. This feature can be attributed to the 
change in the sign of the parameter e across the bounce (similar to the sign change of 5 discussed 

in remark 4). One can write e = ±2H \J~-g^j^ — j^z which suggests that due to the change in the 
sign of mean Hubble rate during the transition from expanding to contracting branch of universe 
along the backward evolution, the parameter e also changes the sign. Owing to this sign change 
in e, although ip does not change, if one starts with an initial condition such that structure near 
the bounce is pancake on one side of the bounce, one does not obtain a pancake structure across. 
In other words, if one starts with e > and if> = 7r/6 forming a pancake like structure near the 
bounce in the expanding branch, the parameters across the bounce will become e < and ip = tt/6 
which is a condition for cigar like structure. 

Fig. [13] shows the behavior of dust filled Bianchi-I spacetime across the bounce. The left plot, 
for which the initial parameter is taken to be e = 34 and ip = tt/6, shows the directional Hubble 
rates. This corresponds to a pancake-cigar transition as discussed above. In the expanding branch, 
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FIG. 13: The left plot shows the directional Hubble rates for a dust filled universe near the bounce when 
the spacetime undergoes a transition from a pancake type structure to a cigar. The right plot shows the 
evolution of the Kasner exponents far from the bounce, all the Kasner exponents tend to take the same 
value (fcj = 2/3) in this limit which is an indication of the isotropization. The initial anisotropy parameters 
for this figure are: e = 34, ip = tt/6 (for the left plot) and ip = 137r/36 (for the right plot). 
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FIG. 14: The left plot shows the directional Hubble rates for a radiation filled universe near the bounce. 
The right plot shows the evolution of Kasner exponents far from the bounce when all the Kasner exponents 
tend to take the same value (i.e. hi = 0.5). This is an indication of isotropization of the spacetime at very 
late time. The initial conditions for this simulation are e = 8.1,?/' = 7r/18. 



two of the directional Hubble rates are zero while the other is positive, this calls for a pancake 
type structure near the bounce. Whereas, on the contracting side across the bounce, two of the 
Hubble rates are negative while the remaining one is positive leading to a cigar like structure. The 
right plot of Fig. [13] depicts the Kasner exponents with the initial anisotropy parameters being 
e = 34 and ip = 137r/36. In this plot it is evident that in the asymptotic limit all the Kasner 
exponents tend to take the same value (i.e. ki —> 2/3) which leads to the isotropization of the 
spacetime such that ctj oc t 2 ^ . This is also a general asymptotic behavior of a classical dust filled 
Bianchi-I, as noted in the previous section. 

An evolution of a radiation Bianchi-I spacetime has been shown in Fig.[T31 Left figure shows 
the directional Hubble rates close to the bounce and the right one portrays the late time behavior 
of the Kasner exponents on the two sides of the bounce. It is evident from the figure that the 
spatial structure near the bounce is cigar on the expanding side as well as on the contracting side 
which entails to a cigar-cigar transition. Similar to the dust model, the pancake-cigar transition 
for radiation also occurs for certain value of ip (ip = tt/6 when e > and ip = tt/2 when e < 0) . The 
asymptotic behavior is evident from the nature of the late time behavior of the Kasner exponents 
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(the right plot in Fig. ll4|) . All the three Kasner exponents tend to take the same value i.e. ki — > 1/2 
which is an indication of the isotropization of the radiation Bianchi-I spacetime in the asymptotic 
limit. This also is a general asymptotic property of a classical radiation filled Bianchi-I universe. 



IV. DISCUSSION 

A key question for various models in LQC is to understand various physical properties in the 
high curvature regime and the way they change in evolution across the bounce. Our goal in this 
manuscript dealt with gaining insights on this issue for Bianchi-I model using effective dynamics. 
Bianchi-I spacetimes are one of the simplest settings where the role of Weyl curvature on the 
nature of the singularities becomes transparent. The dynamics of anisotropic spacetimes has been 
extensively studied in the classical theory In contrast to the isotopic (and homogeneous) 

universe, where the approach to singularity is always point like, a general approach to singularity 
in Bianchi spacetimes can have various geometric structure depending on the anisotropic shear and 
the matter content. These structures, which include barrel, cigar, pancake and a point, are labelled 
using the Kasner exponents in the spacetime metric. Since the formation of these structures is 
tied to the detailed properties of the matter and anisotropics, understanding of their formation 
in LQC, and their evolution across the bounce provides us with valuable hints on the change in 
detailed composition of the universe from one side of the bounce to the other. 

Our analysis is based on the effective spacetime description of LQC, which has been tested 
extensively to confirm with the und erly ing quantu m dy namics for states which lead to a macroscopic 
universe at late times 10 -H, 14, 16, 18, 0, H3, 34]. In the case of Bianchi-I model considered 



here, effective dynamics has provided useful insights on the boundedness of the energ y d ensity, 
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expansion and shear scalars [35j, |43j, |44j and generic resolution of strong singularities [381 ] . The 
modified dynamical equations obtained from the effective Hamiltonian constraint approximate the 
classical dynamical equations to a high order of accuracy at low curvature scales. At these scales, 
two observers starting from the same initial conditions, say evolving backward in the expanding 
branch defined by the mean scale factor, in LQC and in GR measure the same Kasner exponents, 
and thus deduce the same geometric structure which is formed in the high curvature regime. The 
important difference being that in GR structures such as barrel, pancake, cigar and point are 
associated to the big bang singularity, whereas in LQC these structure are finite and non-singular. 
Our detailed analysis for dust, radiation and stiff matter, shows that in the evolution across the 
bounce, the geometric structure in the post-bounce phase can be different from the one in the pre- 
bounce phase. A novel result of our investigation is that there are certain selection rules between 
the transitions of geometric structures across the bounce. These transitions imply change in the 
Kasner exponents across the bounce for asymptotic observers. Many of these Kasner transitions 
occur in a particular range of anisotropic parameters, whereas some occur only for fine tuned values 
(those involving the barrel for stiff matter and the pancake for dust and radiation). Thus, given 
the initial conditions on the anisotropics and the matter, it is possible to identify allowed and 
forbidden transitions. 

We show that there is a hierarchy of different geometrical structures, and in transitions between 
them. Given initial conditions on the matter and anisotropies, we find the existence of allowed and 
forbidden transitions in the effective spacetime description. We find a definite pattern in when a 
particular kind of transition is favored. For example, in the case of stiff matter for large values of 
anisotropic parameter |<5| ~ 1, cigar to cigar transitions are found to be most favored. If the value 
of this parameter (which depends on initial conditions) is chosen to be smaller, the likelihood of 
cigar to cigar transition decreases, and becomes zero for |<5| < 1/3. In a similar way, point like 
to point like transitions are forbidden for \5\ > and are exclusively allowed for |<5| < 1/2. 
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Similar behavior is deduced for all other transitions (summarized in Tables |T] and [TTJ) , and for 
dust and radiation (summarized in Table HV]) . These results clearly show that transitions between 
the anisotropic structures are forbidden in the case of stiff matter, if anisotropic content in the 
universe is small. Similarly, when anisotropies are large, transition between isotropic structures do 
not occur. In summary, depending on the strength of anisotropies, certain transitions are favored 
over others, and some are forbidden. 

Analysis in this article provides useful insights on the structure of the spatial geometry and 
their transitions across the bounce. It brings out some important differences with isotropic models 
in LQC, which arise due to richness of physics introduced by the interplay of the Ricci and Weyl 
curvature components of the spacetime curvature. These Kasner transitions draw a similarity 
with the mixmaster behavior for spatially curved anisotropic models in GR but with some key 
differences. Unlike in GR where numerous Kasner transitions occur in the approach to general 
singularity, in LQC Kasner transition occurs in absence of spatial curvature and only once across 
the non-singular bounce. In an upcoming work, we will provide more details on the BKL like 
phenomena in LQC in a different model [54j]. In future, it will be interesting to analyze the Kasner 
transitions and selection rules using quantum constraint in LQC. Such a study will shed insights on 
the way quantum properties of states affect allowed and forbidden transitions across the bounce. 
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